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ABSTRACT 

Equations  of  motion  are  presented  for  a  three  dimensional  tra¬ 
jectory  simulation  under  the  assumption  that,  during  the  thrust 
phase,  the  configuration  is  a  rigid  body  and  has  90  degree  rota¬ 
tional  symmetry.  For  flight  simulation  during  the  after-burning 
phase  equations  based  on  particle  ballistic  theory  are  presented. 

The  instantaneous  position  of  the  missile  in  space  is  defined 
relative  to  a  spherical,  rotating  earth.  Aerodynamic  coefficients 
considered  during  the  thrust  phase  are  functions  of  Mach  number, 
angle  of  attack  and  effective  roll  angle.  During  the  after-burning 
phase  the  aerodynamic  drag  coefficient  is  a  function  of  Mach  number. 
The  formulation  was  designed  for  the  specific  purpose  of  generating 
accurate  fire  control  data  for  the  ASROC  missile  but  may  have  wide 
application  among  other  types  of  rockets.  For  this  reason  certain 
terms  have  been  included  in  the  model  which  are  not  significant  for 
the  ASROC  application.  Among  them  are  terms  to  simulate  effects  of 
a  rotating  earth.  As  an  option  the  free-flight  after-burning  phase 
may  be  simulated  by  a  six  degree  of  freedom  model. 
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INTRODUCTION 


A  mathematical  model  having  six  degrees  of  freedom,  three  for 
translational  and  three  for  angular  motion,  is  required  for  accu¬ 
rate  simulation  of  the  flight  of  the  ASROC  missile  during  the  thrust 
phase,  while  a  particle  mathematical  model  simulates  with  sufficient 
accuracy  the  motion  after  burning.  Current  theory  regarding  the 
aerodynamic  response  of  ASROC,  as  it  emerges  from  a  closed  cell 
launcher  into  the  air  stream,  specifies  that  the  aerodynamic  moments 
build  up  to  peak  values,  rather  than  instantaneously,  during  the  first 
few  milliseconds  after  launch.  This  effect  has  been  included  in  the 
model. 

At  motor  separation  a  net  forward  thrust  is  exerted  on  the 
second  stage  configuration  producing  disturbances,  including  a 
significant  but  predictable  instantaneous  increase  in  the  velocity 
of  the  second  stage.  Separation  of  the  airframe  from  the  payload 
later  in  the  trajectory  also  causes  disturbances  to  the  motion  of  the 
payload.  In  addition,  if  the  payload  is  a  torpedo  it  is  decelerated 
to  a  safe  water  entry  velocity  by  a  parachute.  (ASROC  trajectory 
events  are  illustrated  in  Figure  1.) 

The  mathematical  model  described  in  this  report  simulates  (l) 
motion  on  a  closed  cell,  non- tip-off  launcher  using  a  one-dimensional 
differential  equation  of  motion  wherein  only  thrust  and  gravity  are 
considered,  (2)  motion  during  the  thrust  phase  using  equations  of 
motion  having  six  degrees  of  freedom  with  the  added  provision  that 
the  components  of  the  restoring  moment  coefficient  generated  by  the 
fins  be  time-dependent  for  a  short  period  after  launch,  and  ( 3) 
motion  during  the  after-burning  phase  using  equations  of  motion  based 
on  particle  ballistics.  Disturbances  accompanying  transition  from 
the  powered  to  the  after  burning  phase  and  from  the  separation  of  the 
airframe  and  the  subsequent  deployment  of  the  parachute  are  considered. 
The  mathematical  model  generates  trajectory  elements  for  each  stage, 
adjusts  the  terminal  position,  velocity  and  time  to  account  for  the 
disturbances,  and  employs  the  adjusted  values  as  initial  conditions 
for  the  succeeding  stage. 

When  the  special  features  incorporated  for  accurate  simulation 
of  the  trajectory  for  the  ASROC  missile  are  removed,  the  model  reduces 
essentially  to  those  described  in  reference  1  as  the  particle  model 
and  the  Littlejohn  model.  The  equations  are  coded,  in  Fortran,  for 
solution  on  either  an  IBM  7090  or  STRETCH  computer  using  th<_ 
Runge-Kutta  procedure  for  numerical  integration. 
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REFERENCE  FRAMES 


a.  Reference  Axes  and  Orientation  Angles 

Two  orthogonal  coordinate  systems  are  used  to  define  the 
missile's  position  and  angular  orientation  for  rigid  body  simula¬ 
tion  during  the  thrust  phase.  The  earth  fixed  reference  system 
(x,y,z)  is  a  right-handed  coordinate  system  with  x  in  a  horizontal 
plane  pointing  down  range;  y  vertical,  positive  upward,  and  z  per¬ 
pendicular  to  the  xy  plane.  The  position  of  the  missile's  center 
of  gravity  is  defined  relative  to  this  system  during  the  powered 
and  the  after-burning  phases. 

The  (X,Y,Z)  reference  system  is  also  an  orthogonal  right- 
handed  coordinate  system  having  its  origin  at  the  center  of  gravity 
of  the  missile.  The  X  axis  is  coincident  with  the  missile's  longi¬ 
tudinal  axis,  positive  forward;  the  Y  axis  is  in  the  vertical  plane 
containing  the  X  axis;  the  Z  axis  is  perpendicular  to  the  XY  plane. 

Orientation  angles  0  and  il i  are  used  to  define  the  instanta- 
eous  elevation  and  azimuth  angles,  respectively,  of  the  missile  with 
respect  to  the  inertial  (x,y,z)  reference  frame  as  shown  in  Figure  2. 
0  is  the  angle  between  the  missile's  longitudinal  axis,  X,  and  the 
horizontal  (xz)  plane  and  \|r  is  the  azimuth  angle  of  the  missile, 
measured  in  the  horizontal  plane,  relative  to  the  x  axis.  0  is  posi¬ 
tive  about  the  positive  Z  axis,  i  is  positive  about  the  negative 
y  axis. 

— ♦ 

The  total  velocity  vector  with  respect  to  the  ground,  V,  is 
computed  from  components  u,  v,  and  w  which  are  components  of  the 
velocity  of  the  origin  of  the  X,Y,Z  axes  relative  to  the  x,y,z 
system,  u,  v,  and  w  are  measured  in  the  direction  of  X,  Y,  and  Z, 
respectively. 

b.  Transformation  of  Reference  Frames 


Velocity  components  of  the  center  of  gravity  u,  v,  and  w 
(along  the  X,Y,Z  axes)  are  expressed  relative  to  the  inertial  (x,y,z) 
reference  system  through  use  of  the  transpose  of  the  transformation 
matrix  A  where 


A  = 


0  cos  i 
9  cos  \lr 


sin  9 
cos  9 
0 


cos  0  sin 
-sin  0  sin 
cos  ii 
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Then 


where  A*  is  the  transpose  of 


the  matrix  A. 


c.  Angular  Velocity  of  the  Reference  Frame 

The  total  angular  velocity  vector  of  the  missile,  w1,  is 
comprised  of  components  p,  q,  and  r  measured  about  the  X,Y,Z  axes, 
respectively,  (see  Figure  2).  Because  the  X,Y,Z  system  does  not 
spin  with  the  missile  (coincident  with  the  constraint  imposed  on 
the  Z  axis  to  lie  in  the  horizontal,  xz,  plane),  the  angular  velo¬ 
city  vector  of  the  missile  reference  frame  is  different  from  the 
angular  velocity  vector  of  the  missile.  It  is  shown  in  Appendix  B 
that  the  angular  velocity  of  the  XYZ  reference  system  presented  in 
the  inertial  reference  system  has  the  components 


(q  tan  6,  -q/cos  6,  r) 


about  the  X,  Y,  and  Z  axes,  respectively.  The  angular  velocity 
vector  of  the  missile  may  then  be  written  as 


where  0  is  the  roll  angle  of  missile  about  the  X  axis  with  respect 
to  the  vertical  plane  containing  the  launching  lugs. 


RIGID  BODY  PHYSICAL  PARAMETERS 


Mass,  moments  of  Inertia  (axial  and  transverse),  and  the  center 
of  gravity  position  are  simulated  as  time  dependent  variables  using 
the  following  equations; 
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a.  Maas 


where  m  Is  the  Instantaneous  mass,  m0  the  Initial  mass,  m  the  rate 
at  which  the  propellant  is  expended,  T  the  Instantaneous  thrust,  and 
V(  the  effective  gas  velocity.  V(  nay  be  expressed  as  the  product 

of  the  specific  Impulse  (lsp)  and  gravity. 

b.  Axial  Moment  of  Inertia 

Jx  *  JXB  +  k*  “p 

where  Ix  is  the  Instantaneous  axial  moment  of  inertia,  IXb  the  axial 
moment  of  inertia  with  all  propellant  expended,  ka  the  square  of  the 
radius  of  gyration  of  the  missile  about  the  longitudinal  axis 

t 

(assumed  to  be  constant),  and  nip  (■  nip0  +•  f  A  dt),  the  Instantaneous 

to 

mass  of  the  propellant. 

c.  Transverse  Moment  of  Inertia 


IY  *  <*2  +  k|  m  -  Xg 

where  Iy  is  the  instantaneous  transverse  moment  of  inertia,  Xg  the 

instantaneous  center  of  gravity  position  with  respect  to  the  missile's 
nose,  k|  the  square  of  the  radius  of  gyration  of  the  propellant  about 

a  transverse  axis  through  the  propellant  center  of  gravity  and  ax 
and  a2  are  constants.  Ihe  method  used  to  determine  the  constants  ax, 
Og,  and  ka  is  given  in  section  d.  below. 

d*  Center  of  Gravity  Position 
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where  Xg  is  the  instantaneous  center  of  gravity  position  with 
respect  to  the  nose;  Xp  is  the  center  of  gravity  position  of  the 
propellant  (assumed  constant). 

The  procedure  for  determining  the  constants  ax,  &2  and 
k|  for  the  ASROC  missile  is  as  follows: 

1.  Obtain  measurements  of  the  center  of  gravity  posi¬ 
tion,  Xg,  for  the  fully  loaded  (unburnt)  round  for  the  round  with 

all  propellant  expended.  Solve  the  system  of  linear  equations 


ai 

xg,o  “  xp  + 


X 


g>b 


+ 


% 


for  ax  and  Xp. 

2.  Obtain  measurements  of  the  transverse  moment  of 
inertia  for  the  fully  loaded  (unbumt)  round  and  for  the  round  with 
all  propellant  expended.  Solve  the  system  of  linear  equations 

!y,o  =  a2  +  kp  mo  -  al  xg,° 

Iy,b  =  a2  +  kf  mb  -  ax  Xg,b 


for  a2  and  k|  . 


EQUATIONS  OF  MOTION 

Components  of  the  total  translational  and  angular  acceleration 
vectors  measured  relative  to  the  X,Y,Z  reference  frame  are,  respec¬ 
tively: 
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where 

Ft  =  force  due  to  thrust 

FA  =  sum  of  the  aerodynamic  forces 

g  =  gravitational  acceleration  vector 
n  x  V  =  acceleration  vector  caused  by  the  earth's  rotation 

i?A  =  sum  of  the  torques  produced  by  aerodynamic  moments 
(including  aerodynamic  trim) 

Mt  *  +  Mj) ,  sum  of  the  torques  produced  by  thrust 

(including  jet  damping  torque,  $j) 

A  differential  equation  representing  the  gravity  free  accelera¬ 
tion  measured  by  an  accelerometer  during  the  thrust  phase  is  also 
computed.  It  is  assumed  that  the  accelerometer  is  oriented  along  the 
X  axis  and  is  therefore  sensitive  only  to  accelerations  along  X.  The 
equation  is 

Vaccel  =  jj  ^  cos  e  +  aSCx(M)  -  (q2  +  r2)  (Xg  -  Xaooel) 

L.  « 

Xaccel  is  the  distance,  along  the  longitudinal  axis  of  the  missile, 
between  the  nose  of  the  missile  and  the  accelerometer. 
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a.  Forces  and  Moments 

Aerodynamic  Forces  (?A) 

Aerodynamic  Forces  considered  are: 

Force  Coefficient 


Drag 

Normal  Force  -  sum  of  the  roll 
and  nonroll  dependent  components 


cx(M,a) 

Cz(M,a)  =  czo(M,a)  + 

CZe(M,a)  sin  4  (A  +  0) 


Magnus 


Cyp(M,a) 


4 


A  is  the  dihedral  angle  between  the  XZ  plane  and  the  yaw  plane  and 
0  is  the  dihedral  angle  between  the  XY  plane  and  a  plane  containing 
the  missile  longitudinal  axis  and  the  launching  lugs.  These  angles 
are  shown  in  Figure  3. 

Mach  number  and  angle  of  attack  are  computed  as  follows: 

Mach  Number: 


M  = 


Angle  of  attack: 


a  =  tan' 


-1  1 


+  w^ 
a  a 


ua 


a. 


tr 


=  tan"1  ^ 


(va  +  ua  tan  5  cos  0tr)2  +  ^wa  +  ua  tan  5  sin  0tr)‘ 


ua 
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where 

Tssl  =  temperature  of  the  air  at  sea  level  under  standard 

atmosphere  conditions 

T*(y)  =  temperature  of  the  air  at  altitude  y 

AgSl  =  speed  of  sound  at  sea  level  under  standard 

atmospheric  conditions 

Va  =  magnitude  of  the  missile  velocity  vector  relative 

to  the  air 

ua,  va,  wa  =  components  of  the  missile's  total  velocity  vector 
relative  to  the  air  measured  along  the  X,Y,Z 
reference  axes,  respectively 

1.  Thrust  Forces  (?T) 


The  mathematical  model  considers  thrust  forces  along  all 
three  missile  fixed  axes,  i.e. ,  forces  resulting  from  a  thrust  vector 
not  coincident  with  the  missile  longitudinal  axis  are  included.  The 
following  examples  are  given  to  illustrate  the  types  of  malalignment 
which  may  be  simulated. 

Example  1.  Linear  Thrust  Malalignment:  The  nozzle  axis 
is  eccentric  with  respect  to  the  missile  X  axis  by  the  distance  rN, 
in  a  direction  defined  by  the  angle  0N. 


Example  2.  Angular  Thrust  Malalignment;  The  axis  normal 
to  the  nozzle  plate  is  malaligned  by  the  angle  e  with  respect  to  the 
missile  X  axis  in  a  plane  oriented  by  the  angle  0O  with  respect  to 
the  body  fixed  Y  axis.  The  axis  of  thrust  passes  through  the  center 
of  gravity  of  the  round. 

Example  3.  Linear  and  Angular  Thrust  Malalignment:  A 
combination  of  examples  1  and  2. 


(The  moments  generated  about  the  body  fixed  axes  due  to  thrust  mal¬ 
alignment  are  discussed  in  the  section  entitled  "Thrust  Moments.") 


The  forces  acting  along  the  missile  fixed  axes  as  a  result  of 
thrust  malalignment  are  given  by  the  equation 


f 
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Vertical  Plane 
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where 


and 


?  =  T 


0, 


th 


1 

0 


^  cos  e 
sin  e  cos  0e 
\  sin  e  sin 


0 

COS  0th 


\  0  sin  0th 


e  =  the  angle  between  the  missile  longitudinal  axis  and 
the  thrust  vector 

T  =  thrust 

0th  *  0th ,o  +  0;  the  instantaneous  dihedral  angle  between 
the  vertical  plane  containing  the  missile  Y  axis  and 
the  plane  containing  the  malaligned  thrust  vector  and 
a  vector  parallel  to  the  missile  longitudinal  axis 

0C  =  the  orientation,  with  respect  to  a  vertical  plane 

formed  by  the  malaligned  thrust  vector  and  the  missile 
longitudinal  axis  of  the  plane  formed  by  the  missile 
longitudinal  axis  and  an  axis  through  the  nozzle  center 

2.  Aerodynamic  Moments  (M^) 


Aerodynamic  moments  considered  are: 

Moment  Coefficient 

Restoring  moment  -  sum  of  the  Cm  =  Cmo(M,a)  +  Cmof  *a)  +  Cmof.  (a) 
body  alone,  booster  fin  and  1  2 

airframe  fin  components  r 

+  +  Cmef2(a)Jcos  4(A+0) 

Side  moment  -  sum  of  booster  r  _  L  ,  c  sin  4  (*  , 

and  airframe  fin  components  n  ~  L  n7*V°  n7f2  a  J  \  r) 
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Magnus  moment  Cnp(M,a) 

Cross  spin  damping  moment  Cmq(M,a) 

Axial  spin  moment  C^gA(M) 

Axial  spin  damping  moment  C^p(M) 


3.  Thrust  Moments  =  Rj  +  #j) 


The  Jet  damping  moment  (ftj)  results  from  an  impingement 
of  gases  inside  the  thrust  chamber  which  exerts  a  damping  moment  on 
the  missile  experiencing  angular  velocities.  The  effect  is  approxi¬ 
mated  by  the  terms 


where  A  is  a  length  associated  with  Jet  damping. 

Thrust  moments  (Ftp)  generated  as  a  result  of  a  mal- 

aligned  thrust  vector  are  also  considered.  Considering  the  cases 
described  in  the  section  on  thrust  forces,  the  moments  corresponding 
to  these  forces  are  given  by  the  vector  equation 

%  -  [fU]  •  [(*  x  f>  +  i?J 

where 

/ -  (L  -  Xg) v 
?  -  j  r„  cos  j 

\rN  sin  / 
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and 


'cos  € 

Mv  =  hT  [  sin  e  cos  0O 
ksin  e  sin  0C 


L  =  length  of  the  rocket 

rN  =  the  perpendicular  distance  between  the  thrust  vector 
and  the  missile  longitudinal  axis 

0N  =  "the  orientation,  with  respect  to  a  vertical  plane  con¬ 
taining  the  missile  longitudinal  axis,  of  the  plane 
containing  rN 


h  =  thrust  axial  moment  coefficient 


kT  *  the  rolling  torque  about  the  thrust  vector  regardless 

of  the  alignment  of  the  latter  and  will  generate  moments 
about  the  three  body  fixed  axes  when  both  e  and  are 
other  than  zero 


The  development  of  terms  to  simulate  a  malaligned  thrust 
vector  are  based  on  the  assumption  of  a  single  nozzle.  However,  the 
equations  apply  to  the  multiple  nozzle  case  provided  that  components 
of  the  thrust  forces  along  the  axes  normal  to  the  nozzle  plate  are 
known.  The  simpliest  case  would  be  one  in  which  the  nozzle  centers 
are  equidistant  from  the  nozzle  plate  center,  have  the  same  cant  angle, 
and  are  spaced  symmetrically  with  respect  to  a  reference  plane  con¬ 
taining  the  nozzle  center. 

4.  Aerodynamic  Trim 

Deviations  in  the  external  configuration  of  the  missile 
from  the  symmetrical  shape  to  which  basic  wind  tunnel  aerodynamic  data 
apply  may  result  in  the  aerodynamic  normal  force  and  restoring  moment 
being  zero  for  a  nonzero  angle  of  attack.  The  magnitude  of  this  angle 
of  attack  is  referred  to  as  the  trim  angle,  6.  The  axis  defined  by  S 
is  the  aerodynamic  trim  axis  which  specifies  the  direction  of  the  air¬ 
speed  vector  for  which  there  are  no  normal  forces.  The  resulting 
cross  velocity  components  with  respect  to  the  air  under  conditions  of 
aerodynamic  trim  are 
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5^0 


+  [  ua  tan  5  cos  $tr 

ua  tan  5  sin  $tr 


The  orientation  angle  0tr  is  the  instantaneous  dihedral  angle  about 

the  missile  X  axis  between  the  XY  plane  and  the  plane  containing  the 
aerodynamic  trim  axis  and  the  missile  X  axis.  At  any  instant 

0tr  =  0tr,o  +  0 

where  the  subscript  "0"  indicates  initial  value. 


b.  Gravity 


The  simulation  of  acceleration  due  to  gravity  is  based  on 
the  inverse  square  lav  of  force  as  applied  to  a  spherical  earth: 


?  =  - 


So 


r 


3 


where 


gQ  =  the  gravitational  acceleration  at  the  origin  of  the 
earth  fixed  reference  system 

r0  =  the  radius  of  the  earth 

r  =  the  position  vector  of  the  missile  center  of  gravity 
in  space 

r  =  |?| 


At  the  origin  of  the  earth  fixed  (x,y,z)  system  the  gravity  vector 
is 


At  any  point  on  the  trajectory,  away  from  the  origin,  the  gravity 
vector  is  approximated  by 

2?  =  <?  +  AS? 

6approx 
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where 

£&0  =  the  derivative  of  ?  (with  respect  to  position) 

evaluated  at  the  origin  multiplied  by  the  change 
in  position  (ag0  =  vg0  Ar) 

For  the  present  case,  the  z  component  of  position  is  small  enough 
to  be  neglected  insofar  as  its  influence  on  £  is  concerned.  In  the 
simulation,  therefore,  the  following  gravity  vector  is  employed: 


Note  that  elsewhere  in  the  report  the  radius  of  the  earth  is  denoted 
by  re,  r  is  the  component  of  the  total  angular  velocity  measured  about 

the  Z  axis  and  the  subscript  "0"  denotes  time  of  ignition  of  the  pro¬ 
pellant  . 

c.  Coriolis  Acceleration 

The  acceleration  vector  resulting  from  the  assumption  of  a 
rotating  earth  presented  in  the  (x,y,z)  reference  system,  is 


where  n  is  the  singular  velocity  vector  of  the  earth  (directed  along 
the  polar  axis)  and  the  total  velocity  vector  of  the  missile 
relative  to  the  earth  fixed  system.  Expressions  for  the  components 
of  J?  are  derived  as  follows: 

Define  £  as  the  angle  between  the  earth's  axis  and  the 
horizontal  plane  (positive  in  the  northern  hemisphere)  and  AzL  as 

the  angle  from  north  to  the  horizontal  projection  of  the 
launcher  or  the  azimuth  angle  of  the  positive  x  axis.  The 
positive  z  axis  will  then  have  azimuth  90°  +  AzL.  Therefore, 
a  unit  vector  drawn  horizontally  northward  will  have  x  and  z 
components  which  are,  respectively  cos  AzL  and  cos 
(90°  +  Azl)  *  -  sin  Azl.  The  horizontal  projection  of  j?  points 
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northward  and  has  length  fi  cgs  L,  so  its  x  and  z  components  are 
fl  cos  L  cos  Azl  and  -  fi  cos  L  sin  AzL ,  respectively.  Its  y  com¬ 
ponent  is  o  sin  L.  The  three  components  of  J ?  may  then  be  written 
as 


Then, 


/  cos  L  cos  Azl^ 

/  ax  \ 

-Si 

fit  = 

sin  L 

- 

®y 

-cos  L  sin  Az. 

a* 

\  / 

\  1 

Presentation  of  the  vector  in  the  missile  fixed  reference  system 
is  accomplished  by  premultiplying  the  above  expression  t>y  the 
A  matrix . 


d.  Atmosphere 
Wind 


Horizontal  wind  is  expressed  as  a  velocity,  W,  and  an 
orientation,  i|rw,  where 


=  Azw  -  Azl 


Azw  is  the  angle  between  north  and  the  direction  from  which  the 

wind  blows.  AzL  is  the  angle  between  north  and  the  vertical  plan-. 

through  the  launcher.  Thus,  \Jrw  is  the  angle  between  the  vertical 

plane  containing  the  launcher  line  and  the  wind  vector.  Components 
of  W  expressed  relative  to  the  (x,y,z)  axes  are  then  given  by 


,  w_\  / cos  (Azw  -  Azl) 

l  \  j 


W, 


*  -  w 


^  sin  (Azw  -  Azl) 
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The  computer  program  provides  for  the  computation  of  the  wind 
vector  from  components  of  ambient  wind  and  from  a  moving  launcher 
as  follows: 


W» 


/ 


\ 


cos  (Azw^ 
0 

sin  (Azw 


(Azt 

0 

(  Azt 


where 


WA  =  magnitude  of  the  ambient  wind 


Azw  *  angle  between  north  and  the  direction  from  which  the 
A  wind  blows 

W0  =  velocity  of  the  launcher  (or  a  moving  cart  on  which 
the  launcher  is  mounted) 

Azt  «  angle  between  north  and  the  direction  of  motion  of  the 
launcher  (or  cart) 


Since  the  aerodynamic  forces  and  moments  on  the  missile  are  functions 
of  missile  velocities  relative  to  the  air,  the  components  (ua,  va,  wa) 
measured  along  the  missile  reference  axes  are  computed  as  follows: 


It  follows  that  the  total  velocity  relative  to  the  air  mass,  Va,  is 

Va  -  (uf  +  vf  +  vf)1/2 

The  air  density,  p,  and  air  temperature,  T*,  are  inputs  as  functions 
of  altitude,  y,  and  must  be  provided  in  the  special  format  and  size 
described  in  Appendix  C. 

Force,  moment  and  acceleration  terms  contributing  to  the  total 
translational  and  angular  acceleration  vectors  are  summarized  in 
Table  1. 
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INITIAL  CONDITIONS  AND  STAGING  CONSIDERATIONS 

Hie  initial  conditions  for  the  missile  trajectory  at  launch 
presented  in  the  Inertial  reference  system  are: 

a.  For  the  position  vector,  !? 


where  R  is  the  length  and  Q  the  elevation  angle  of  the  launcher, 
b.  For  the  velocity  vector,  S 


/  x  \  /  cos  e  \ 

y  ) =  R  l  sin  0  I 

\J  No  ' 

where  the  vector  f?  is  obtained  by  integrating  the  one  dimensional 
differential  equation 


g  sin  6 


to  the  time  when  R  equals  the  launcher  length.  The  differential 
equation  used  for  simulating  the  gravity  free  acceleration  sensed 
by  the  accelerometer  during  launch  is 

V  i  ®  — 
vaeoel  -  £ 


FREE  FLIGHT  PARTICLE  TRAJECTORY  MODEL 


Impact  predictions  based  on  a  simulation  of  missile  motion 
using  a  three  dimensional  particle  model  between  thrust  termina¬ 
tion  and  impact  were  almost  as  accurate  as  would  have  been  obtained 
using  rigid  body  equations  of  motion.  Hie  above  consideration  and 
the  resulting  reduction  in  computer  time  make  the  particle  model 
a  logical  choice  for  simulating  the  after  burning  phase. 
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In  the  particle  model  the  total  translational  acceleration 
vector,  expressed  in  component  form  relative  to  the  inertial 
(x,y, z)  reference  system,  is  given  by 


To  start  computations  of  the  second  stage  of  the  trajectory, 
initial  values  of  time,  position  and  velocity,  appropriate  for  the 
time  at  which  maximum  velocity  is  attained,  must  be  generated  from 
the  terminal  point  computed  using  the  six  degree  of  freedom  model. 
The  procedure  coded  to  obtain  the  initial  values  is  as  follows: 

a.  Compute  the  time  maximum  velocity  is  expected  to  occur 
by  adjusting  the  time  of  motor  separation  by  a  correction,  based 
on  observed  performance,  supplied  as  input, 

b.  compute  components  of  velocity  appropriate  at  the  time 
maximum  velocity  is  expected  to  occur  by  adjusting  values  computed 
using  the  six  degree  of  freedom  model  at  time  of  motor  separation 
by  amounts,  based  on  observed  performance,  supplied  as  inputs,  and, 

c.  compute  position  coordinates  appropriate  for  the  time 
maximum  velocity  is  expected  to  occur  by  adjusting  the  coordinates 
computed  using  the  six  degree  of  freedom  model  at  motor  separation 
by  an  amount  approximated  using  the  assumption  of  gravity  free 
travel  in  a  vacuum. 

The  equations  coded  to  perform  these  functions  are: 

a.  Time 


ty  =  t«  o  i  ~ 

where  is  the  tljne  maximum  velocity  is  attained  and  tco  is  the 

motor  separation  time.  At1-2  is  an  input  constant  and  represents 

the  average  time  between  separation  of  the  motor  from  the  payload 
and  the  time  of  maximum  velocity  of  the  payload. 
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b.  Velocity 


where  Ax  and  Ay  are  the  x  and  y  components  of  velocity  resulting 
from  the  impulsive  force  Imparted  to  the  second  stage  at  motor 
separation.  Ax  and  Ay  are  evaluated  as  functions  of  the  expected 
velocity  at  motor  separation  and  supplied  as  Inputs.  Mass  and 
second  stage  diameter  are  the  only  other  new  Inputs  required. 

c .  Position  Coordinates 


Motion  between  the  time  that  separation  of  the  airframe  is 
Initiated  and  the  time  when  the  parachute  is  fully  deployed  is  not 
simulated.  The  initial  values  of  time,  velocity,  and  position 
necessary  to  "restart”  the  trajectory  at  time  of  full  parachute 
deployment  are  obtained  by  adjusting  values  generated  at  time  of 
airframe  separation  as  follows: 

a.  Compute  estimate  o  f  time  required  for  the  parachute  to 
became  fully  deployed,  and  adjust  the  time  of  airframe  separation 
by  this  estimate. 

b.  Compute  position  coordinates  appropriate  for  adjusted  time 
by  applying  estimates,  based  on  the  assumption  of  gravity  free  travel 
between  the  time  of  airframe  separation  and  the  adjusted  time,  com* 
puted  as  in  a.  above.  It  is  assumed  that  there  is  no  loss  in  veloc¬ 
ity  during  parachute  deployment,  (if  the  configuration  is  a  depth 
charge,  the  parachute  is  not  used.  It  is  assumed  that  separation  of 
the  airframe  is  Instantaneous  and  that  no  disturbances  are  imparted 
to  the  payload.) 
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The  equations  coded  to  perform  these  functions  are: 
a.  Time 


3,0 


=  t  +  At 


2-3 


where 


t3j0  =  time  when  the  parachute  is  fully  deployed. 

t2  =  time  that  separation  of  the  airframe  is  initiated. 

At2_3  =  Cx  +  C2  Vt2  +  C3(Vt2)2  ,  scaled  value,  based  on 

observed  performance,  of  the  time  required  for 
full  deployment  of  the  parachute.  (Cx,  C2  and 
C3  *  0  when  flight  of  the  depth  charge  is 
simulated . ) 

Vt  =  missile  velocity  at  time  of  airframe  separation, 
b .  Position  Coordinates 


Complete  descriptions  of  separation  phenomena  and  disturbances 
during  separation  of  the  airframe  and  the  subsequent  deployment  of 
the  parachute  are  given  in  reference  2. 
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GLOSSARY 


(A)  transformation  matrix  relating  the  missile  axis 

reference  system  to  the  earth-fixed  reference 
system 


Azt 


standard  sea  level  value  for  the  speed  of  sound 

azimuth  angle  of  the  launcher  (measured  clock¬ 
wise  from  north) 

azimuth  angle  of  the  track  (measured  clockwise 
from  north  to  the  direction  of  motion  of  the  cart) 


azimuth  angle  of  the  ambient  wind  (measured  clock¬ 
wise  from  north  to  the  direction  from  which  the 
wind  blows) 


(ax,  ay,  az)  components  of  the  earth's  angular  velocity  pre¬ 
sented  in  the  earth-fixed  reference  system 


(C-l,  C2,  C3)  coefficients  for  the  quadratic  equation  used  to 
compute  parachute  deployment  delay  time 


(^i*  ,  d3) 


diameters  of  the  first,  second,  and  third  stage 
configurations,  respectively 


F 


total  external  force  exerted  on  the  rocket 


(f*,  f2) 


as  subscripts  denote,  respectively,  booster  and 
airframe  fin  contributions 


components  of  !?[  along  the  missile  fixed  X,Y,Z 
axes,  respectively 


g  acceleration  due  to  gravity 

Ix  instantaneous  axial  moment  of  inertia  of  the 

rocket 


instantaneous  transverse  moment  of  inertia  of 
the  rocket 


k2 


squahe  of  the  radius  of  gyration  of  the  rocket 
about  the  longitudinal  axis  of  the  rocket 
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GLOSSARY  (Continued) 

square  of  the  radius  of  gyration  of  the  propellant 
about  a  transverse  axis  through  the  propellant 
center  of  gravity 

l  length  associated  with  jet  damping  moment 

L  length  of  the  missile 

L  latitude  of  the  launcher 


m 


0 

(P>  r) 


R 


r 


e 


r 


N 


s 


instantaneous  mass  of  the  rocket  and  propellant 

instantaneous  mass  of  the  propellant 

as  a  subscript  denotes  time  of  propellant  ignition 

components  of  the  instantaneous  angular  velocity  of 
the  missile  with  respect  to  the  inertial  reference 
axes,  measured  about  the  missile  fixed  (X,Y,Z)  axes, 
respectively 

length  of  launcher 

radius  of  the  earth 

perpendicular  distance  between  the  malaligned  thrust 
vector  and  the  missile  center  of  gravity  position 

maximum  cross  sectional  area  of  the  missile 


t 


'max 


At 


1-2 


time  from  ignition 

time  maximum  velocity  is  reached 

time  of  full  deployment  of  the  parachute 

time  between  motor  separation  and  time  of  maximum 
velocity 


^2-3 


time  between  separation  of  the  airframe  from  the 
payload  and  full  deployment  of  the  parachute 


T 


thrust 
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T*(y)  temperature  of  the  air,  function  of  altitude,  (y) 


(u,  v,  w)  components  of  the  instantaneous  translational 

velocity  of  the  missile  center  of  gravity  with 
respect  to  the  earth,  presented  along  the  X,Y,Z 
axes,  respectively 

( ua ,  va,  va)  components  of  the  instantaneous  translational 
velocity  of  the  missile  center  of  gravity  with 
respect  to  the  air,  presented  along  the  X,Y,Z  axes, 
respectively 

V  magnitude  of  the  instantaneous  translational  veloc¬ 

ity  vector  of  the  missile  with  respect  to  the  earth 

Va  magnitude  of  the  instantaneous  velocity  of  the 

missile  with  respect  to  the  air 

Vaccel  integral,  with  respect  to  time,  of  the  gravity  free 

acceleration  of  the  missile  fixed  accelerometer 

Vg  average  velocity  of  the  jet  gases 

WA  magnitude  of  the  velocity  vector  of  the  ambient  wind 

with  respect  to  the  earth 

Wc  cart  speed,  relative  to  the  earth 

(Wx,  WY,  W2)  components  of  the  wind  vector  resolved  along  the 
X,Y,Z  axes,  respectively 

(Wx,  Wy,  W2)  components  of  the  wind  vector  resolved  along  the 
x,y,z  axes,  respectively 

(x,y,z)  earth-fixed  reference  frame  (an  orthogonal,  right- 

handed  system  of  axes  having  x  horizontal  and  point¬ 
ing  down  range,  y  vertically  upwards,  and  z  per¬ 
pendicular  to  the  xy  plane) 

(X,Y,Z)  missile  axis  reference  system  (an  orthogonal,  right- 

handed  system  of  axes  with  X  coinciding  with  the 
longitudinal  axis  of  the  rocket;  Y  is  in  the  vertical 
plane  and  perpendicular  to  X;  Z  is  constrained  to  lie 
in  the  horizontal  (xz)  plane 


3 


a 

ax,  as 

8 


8a 

€ 


input  constants  used  to  generate  the  instantaneous 
center  of  gravity  position  and  the  instantaneous 
transverse  moment  of  inertia  of  the  missile 


angle  of  trim  (angle  between  the  missile  X  axis 
and  the  aerodynamic  trim  axis,  the  latter  being 
a  missile  fixed  axis  which  defines  the  direction 
of  the  airspeed  vector  for  which  there  are  no  aero¬ 
dynamic  normal  forces) 

fin  cant  angle  (dihedral  angle  between  the  plane 
of  the  fin  and  the  plane  containing  the  longitudinal 
axis  of  the  missile  and  the  center  of  the  fin) 


angle  between  the  missile  X  axis  and  the  thrust 
vector 


6  angle  of  the  missile  longitudinal  axis  with  respect 

to  the  horizontal  (xz)  plane  (positive  toward  the 
positive  y  axis) 

K  thrust  axial  moment  coefficient 

Xg  distance  between  the  nose  and  the  instantaneous 

center  of  gravity  position  of  the  missile 

Xaocei  distance  between  the  nose  of  the  missile  and  the 

center  of  gravity  of  the  accelerometer 

Xp  distance  between  the  nose  of  the  missile  arid  the 

center  of  gravity  of  the  propellant 

X  (a)  distance  between  the  nose  of  the  missile  and  the 

"  '  center  of  pressure  of  the  normal  force  for  the 

"body  minus  fins"  configuration  of  the  missile 

XMi,  XN2  distance  between  the  nose  of  the  missile  and  the 

centers  of  pressure  of  the  normal  forces  generated 
on  the  booster  and  airframe  fins,  respectively 

A  the  dihedral  angle  between  the  XZ  plane  and  the  yaw 

plane 


i 


i 
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GLOSSARY  (Continued) 

Hx(t),  |i2(t)  time  dependent  scaling  factors  for  side  and  restoring 
moment  coefficients  generated  by  the  booster  and 
airframe  fins,  respectively 

density  of  the  air  (function  of  altitude,  y) 

j } i  missile  roll  angle,  measured  about  the  missile  X  axis 

0C  dihedral  angle  between  the  plane  formed  by  the  thrust 

vector  and  the  missile  axis  and  a  vertical  plane, 
parallel  to  the  vertical  plane  containing  the  missile 
longitudinal  axis,  through  the  center  of  the  nozzle 
base 

0N  dihedral  angle  between  the  vertical  plane  containing 

the  missile  longitudinal  axis  and  the  plane  formed  by 
the  missile  longitudinal  axis  and  the  malaligned 
thrust  vector 

0th  instantaneous  dihedral  angle  between  the  missile 

Y  axis  and  the  plane  containing  the  malaligned  thrust 
vector  and  a  vector  parallel  to  the  missile  X  axis; 
positive  when  measured  clockwise  from  the  vertical 
plane 

0tr  instantaneous  dihedral  angle  between  the  missile  Y 

axis  and  the  plane  containing  the  aerodynamic  trim 
axis  and  a  vector  parallel  to  the  missile  X  axis; 
positive  when  measured  clockwise  from  the  vertical 
plane 

\j r  azimuth  angle  of  the  missile  measured  in  the 

horizontal  plane  relative  to  the  x  axis  (positive 
about  the  negative  y  axis) 

\|fw  angle  between  the  horizontal  wind  component  and  the 

x  axis  (positive  about  the  negative  y  axis) 

u  total  angular  velocity  of  the  missile  with  respect 

to  the  earth 

q  angular  velocity  of  the  earth  about  the  polar  axis 
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CiP(M) 

WM> 

Cjj|,(M,a) 

Cnp(M>Ql) 

CN0f1(O!) 

CNOf2^a^ 

CN8f2^a^ 

CY7f1^Q!^ 


GLOSSARY  (Continued) 

spin  damping  moment  coefficient  (function  of  Mach 
number) 

spin  moment  coefficient  due  to  cant  angle  of  the 
fins  (function  of  Mach  number) 

Normal  force  coefficient  for  the  configuration 
without  fins  (function  of  Mach  number  angle  of 
attack) 

magnus  moment  coefficient  (function  of  Mach  number 
and  angle  of  attack) 

nonroll-dependent  component  of  the  component  of  the 
total  normal  force  on  the  booster  fins  (function 
of  yaw) 

nonroll -dependent  component  of  the  component  of  the 
total  normal  force  coefficient  on  the  airframe  fins 
(function  of  yaw) 

roll-dependent  component  of  the  component  of  the 
total  normal  force  on  the  booster  fins  (function  of 
yaw) 

roll  dependent  component  of  the  component  of  the 
total  normal  force  on  the  airframe  fins  (function 
of  yaw) 

side  force  coefficient,  booster  fin  component 
(function  of  yaw) 


side  force  coefficient,  airframe  fin  component 
(function  of  yaw) 


Cnq(M,a)  cross  spin  damping  moment  coefficient  (function  of 

Mach  number  and  yaw) 

Cx(M,a)  axial  drag  force  coefficient  for  the  first  stage 

configuration  (function  of  Mach  number  and  yaw) 

C  (M)  axial  drag  force  coefficient  for  the  second  stage 

configuration  (function  of  Mach  number) 
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Cyp(M,a) 

Cz0(M,a) 


cza(M,a) 


axial  drag  coefficient  for  the  third  stage  con¬ 
figuration  (function  of  Mach  number) 


magnus  force  coefficient  (function  of  Mach  number 
and  angle  of  attack) 


nonroll-dependent  component  of  the  total  normal 
force  coefficient  (function  of  Mach  number  and 
yaw) 


roll-dependent  component  of  the  total  normal  force 
coefficient  (function  of  Mach  number  and  yaw) 
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DERIVATION  OF  ANGULAR  VELOCITY  COMPONENTS  OF 
THE  MISSILE  RELATED  TO  THE  INERTIAL  REFERENCE  FRAME 

Since  the  Z  axis  of  the  missile-fixed  reference  frame  is  con¬ 
strained  to  lie  in  a  horizontal  plane,  consideration  must  be  given 
to  the  relationship  between  the  angular  velocity  components  of  the 
missile  and  the  angular  velocity  of  the  coordinate  system.  Refer¬ 
ring  to  the  figure  below,  if  we  let  e,,  e_,  and  e_  be  unit  vectors 

along  the  X,  Y,  and  Z  axes,  respectively,  the  angular  velocity  of 
the  e^  e2,  e3  system  will  have  the  component  6  about  the  e3  and 

-  i'  about  the  upward  vertical,  e2.  Let  Q  =  angular  velocity  of 

ex,  e2,  e3  with  respect  to  the  x,y,z  axes.  Then 

n  =  Q  e3  -  if  e2 

Resolved  in  the  instantaneous  directions  of  the  unit  vectors 


Now,  let  w  be  the  total  angular  velocity  of  the  rocket  such  that,  in 
component  form 


Since  ex,  e2,  e3  do  not  turn  about  ex,  its  angular  velocity  is 
independent  of  the  p  component.  Because  e3  is  always  in  a  horizontal 
plane  q.  and  are  the  same  about  the  y  and  e2  axes,  respectively,  as 
are  the  angular  velocities  r  and  n  about  the  e3  axis.  That  is, 

=  q 

n3  =  r 


l 
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Further 

qx  ■  -q  sin  0 

-  (-q/cos  6)  sin  0 

\ 

■  -q  tan  0  . 

Therefore, 


/ 

(  \ 

/  * 

fp  -  q  tan  0 

• 

* 

m 

-q/cos  0 

w 

r  / 

Note  that  elsewhere  in  the  report  n  represents  the  angular  rate  of 
rotation  of  the  earth  about  its  polar  axis. 
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INPUT  PARAMETERS  AND  INITIAL  CONDITIONS  REQUIRED  FOR 
TRAJECTORY  COMPUTATIONS  ON  THE  IBM  7090  COMPUTER 


Exponential  Input  (Format  5E  14.7) 


t0  (sec) 
m0  ( slugs) 
-R  (ft) 
g  (ft/sec2) 
0  (deg) 

Azw  (deg) 


WA  (ft/sec) 


&uL 

Avl  (ft/sec) 
Avl 


> 


PL 


(rad./ sec) 


Vg  (ft/ sec) 
constant 


motor  ignition  time 

mass  at  motor  ignition  time 

negative  of  the  launcher  length 

gravitational  acceleration 

launcher  elevation  angle 

azimuth  angle,  with  respect  to 
north,  of  the  wind  (measured  from 
north  to  the  direction  from  which 
the  wind  comes) 

wind  velocity  when  constant  winds 
for  all  altitudes  are  assumed 


incremental  velocity  components 
introduced  at  time  of  launch 


missile  angular  velocity  components 
at  launch 


effective  gas  velocity 
scaling  factor  for  thrust 


t  (sec)  or  V  (ft/sec)  value  sensed  to  stop  computations  in 

the  first  stage 

constant  scaling  factor  for  air  density 

0L  (deg)  roll  angle  of  the  missile  fins 
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(deg) 


azimuth  angle  of  the  missile  at 
launch 


yo  (ft)  altitude  at  ignition 

control  word  for  wind  if  zero,  the  first  wind  table  is 

used;  if  nonzero,  the  next  wind 
table  is  used 


control  word  for  air  tempera-  if  zero,  the  first  air  temperature 
ture  and  air  density  and  air  density  table  is  used;  if 

nonzero,  the  next  air  temperature 
and  air  density  table  is  used 


Azt  (deg) 
Azl  (deg) 


azimuth  angle  of  the  track  (measured 
from  north  to  the  direction  of  motion 
of  the  cart) 

azimuth  angle  of  the  launcher  with 
respect  to  north 


1  (ft) 


length  associated  with  Jet  damping 


fi  (rad/sec) 


angular  rate  of  rotation  of  the  earth 
about  its  polar  axis 


IXB  (slug- ft2) 


missile  axial  moment  of  inertia  when 
all  propellant  is  expended 


t  (sec) 


L  (ft) 
Xp  (  ft ) 


^■accel  (f^J 


time  when  the  scaling  factors  for  the 
restoring  moment  coefficient,  |ix  2(t), 
equal  unity 

integration  interval  used  between 
time  where  u.  _(t)  =1  and  the  end  of 

the  first  stage 

length  of  the  missile 

distance  from  the  nose  of  the  missile 
to  the  center  of  gravity  of  the  pro¬ 
pellant 

distance  from  the  nose  of  the  missile 
to  the  center  of  gravity  of  the 
accelerometer 
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k2  (ft2) 
l/r.  (ft'1) 

L  (deg) 

Wc  (ft/sec) 

€  (deg) 

0N  (deg) 


0th, o  (deg) 
0tr,o  (deg) 

5  (deg) 
k  (ft)  • 

C*!  (slug-ft) 

a2  (slug-ft2) 

IDpo  (Slugs) 
AtL  (sec) 

Atx  (sec) 

8a  (rad) 
d2  (ft) 


radius  of  gyration  squared 

reciprocal  of  the  radius  of  the 
earth 

latitude  of  the  launcher 

cart  velocity  (used  when  simu¬ 
lating  firings  from  a  moving  cart) 

angle  between  the  mal aligned  thrust 
vector  and  the  missile  longitudinal 
axis 

orientation  with  respect  to  the 
vertical  plane  containing  the  launch¬ 
ing  lugs  of  the  plane  containing  rN 

orientation  of  the  plane  containing  e 

orientation  of  the  plane  containing  5 

aerodynamic  trim  angle 

length  associated  with  thrust  axial 
moment 

constant  used  to  define  the  instanta¬ 
neous  center  of  gravity  position  and 
the  transverse  moment  of  inertia 

constant  used  to  define  the  instanta¬ 
neous  transverse  moment  of  inertia 

mass  of  the  propellant  at  ignition 

integration  interval  to  be  used  when 
simulating  motion  on  the  launcher 

integration  interval  to  be  used  when 
simulating  motion  during  first  stage 

cant  angle  of  the  fins 

diameter  of  the  first  stage 
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AT» 


addend  used  to  convert  temperature 
from  unit  used  in  the  table  to  °C 


linear  unit  used  per  foot  conversion  factor  to  convert  linear 

and  velocity  units  given  in  the  input 
to  ft  and  ft/sec 

angular  unit  used  per  radian  conversion  factor  to  convert  angular 

units  given  in  input  to  rad  and 
rad/sec 


angular  unit  used  per  degree 

k2  (ft2) 

ta  (sec) 

constant 
r„  (ft) 

0o  (deg) 


constant 

constant 

constant 

constant 

constant 

constant 

constant 


conversion  factor  to  convert  angular 
unit  used  in  computations  to  deg  and 
deg/ sec 

square  of  the  radius  of  gyration  of 
the  propellant 

time  at  which  the  accelerometer  starts 
performing 


scaling  factor  for  Cyp 

perpendicular  distance  between  the 
mal aligned  thrust  vector  and  the 
missile  center  of  gravity  position 

dihedral  angle  between  the  plane  formed 
by  the  thrust  vector  and  the  missile 
axis  and  a  vertical  plane,  parallel  to 
the  vertical  plane  containing  the 
missile  longitudinal  axis,  through  the 
center  of  the  nozzle  base 

scaling  factor  for  CNofi 

scaling  factor  for  Cj|0f2 

scaling  factor  for  CNafi 

scaling  factor  for  Cjjef ^ 

scaling  factor  for  C _ 

YTf  1 

scaling  factor  for  CY7f£ 
scaling  factor  for  Cmq 
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constant 

constant 

constant 

constant 

constant 

constant 

constant 

constant 

V  (ft> 

1^  (slug-ft2) 
1^  (slug-ft2) 
m2  (slugs) 
d2  (ft) 
t2E  (sec) 

Ax2  (ft/sec) 

A/2  (ft/sec) 

Azg  (ft/sec) 


scaling  factor 

scaling  factor 

scaling  factor 

scaling  factor 

scaling  factor 

scaling  factor 

scaling  factor 

scaling  factor 

distance  from  the  nose  to  the 
center  of  gravity  of  the  second 
stage  missile  configuration 

axial  moment  of  inertia  of  the 
second  stage  configuration 

transverse  moment  of  inertia  for 
the  second  stage  configuration 

mass  of  the  second  stage  configura¬ 
tion 

diameter  of  the  second  stage  con¬ 
figuration 

time  vhen  separation  of  the  air¬ 
frame  is  initiated 

horizontal  component  of  impulsive 
velocity  imparted  to  the  second 
stage 

vertical  component  of  impulsive 
velocity  imparted  to  the  second 
stage 

lateral  component  of  impulsive 
velocity  imparted  to  the  second 
stage 


for  ^(t) 
for  n2(t) 
for  Cz0 

for  Cze 
for  CNb 

for  Cip 

f°r 

for  Cv 
Ai 
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Atlg  (sec) 


t2  (sec) 


x2  (ft)  ■'j 

y2  (ft)  l 

z2  (ft)  J 


x2  (ft/sec) 
y2  (ft/sec) 
z2  (ft/sec) 


*2X1 

^£2 


Ata  ( sec) 


constant 

m3  (slugs) 
d3  (ft) 


time  between  motor  separation  and 
the  time  maximum  missile  velocity 
is  attained 

time  maximum  velocity  is  attained 
(used  only  when  starting  a  tra¬ 
jectory  in  the  second  stage) 


position  coordinates  of  the  second 
stage  (used  only  when  starting  a 
trajectory  in  the  second  stage) 


components  of  velocity  for  the  second 
stage  (used  only  when  starting  a  tra¬ 
jectory  in  the  second  stage) 


bracketing  times  used  to  determine 
the  time  of  parachute  deployment 
which  would  minimize  third  stage 
flight  time  while  satisfying  a  crit¬ 
ical  water  entry  velocity 

integration  interval  to  be  used  when 
simulating  motion  during  the  second 
stage  (particle  model) 

scaling  factor  for  the  second  stage 
drag  coefficient 

mass  of  the  third  stage  configuration 

diameter  of  the  third  stage  configura¬ 
tion 


y  (ft) 


descending  altitude  on  which  to  stop 
the  computations 


constant 


scaling  factor  for  the  third  stage 
drag  coefficient,  CX3(m) 


6 


NWL  REPORT  NO.  1812 


c2 

c3 


(sec) 

( sec2/ft) 
( sec3/ft2 


t3  ( sec) 


x3  (ft) 

y3  (ft) 
z3  (ft) 


j 


y3 

*a 


(ft/sec) 

(ft/sec) 

(ft/sec) 


At3  (sec) 


VWE  (ft/sec) 


AV 


WE 


constants  for  the  quadratic  equa¬ 
tion  used  to  evaluate  At„  the 

2-3' 

parachute  delay  deployment  time 


time  when  the  parachute  is  fully 
deployed  (used  only  when  starting 
a  trajectory  in  the  third  stage) 


position  coordinates  when  the  para¬ 
chute  is  fully  deployed  (used  only 
when  starting  a  trajectory  in  the 
third  stage) 


components  of  velocity  when  the 
parachute  is  fully  deployed  (used 
only  when  starting  a  trajectory  in 
the  third  stage) 


integration  interval  to  be  used  when 
simulating  motion  during  the  third 
stage 

water  entry  velocity  (use  zero  except 
when  running  the  special  program  t  o 
minimize  third  stage  flight  time  with 
the  parachute  deployed) 

tolerance  on  the  water  entry  velocity 
above 


6  (deg)  elevation  angle  of  the  missile  at  the 

2  beginning  of  the  second  stage 

i) f  (deg)  azimuth  angle  of  the  missile  at  the 

beginning  of  the  second  stage 
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Ap2  (rad/sec) 
dcig  (rad/sec) 
dr2  (rad/sec) 


Incremental  components  of  missile 
angular  velocity  for  the  second  stage 
(used  only  when  starting  computations 
in  the  second  stage  using  rigid  body 
equations) 


p2  (rad/sec) 
q2  (rad/sec) 
r2  (rad/ sec) 


components  of  missile  angular  velo¬ 
city  for  the  second  stage 


6Aa  (rad)  cant  angle  of  the  fins  of  the  second 

stage  configuration 


constant 

scaling 

factor 

for 

cx2 

constant 

scaling 

factor 

for 

CZ0,2 

constant 

scaling 

factor 

for 

CZ8,2 

constant 

scaling 

factor 

for 

c 

me, 2 

constant 

scaling 

factor 

for 

CNp  1 2 

constant 

scaling 

factor 

for 

^m,2 

constant 

scaling 

factor 

for 

n 

^mq,2 

constant 

scaling 

factor 

for 

cyp>2 

constant 

scaling 

factor 

for 

C-lp,2 

constant 

scaling 

factor 

for 

c«.,2 

constant 

scaling 

factor 

for 

Cy7>2 

yaE  altitude  (descending)  on  which  the 

second  stage  is  terminated 
( option) 
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Integral  Input  (Format  1013) 


Locations 
1,  2,  3 


4,  5,  6 


7,  8,  9 
10,  11,  12 


13,  14,  15 


Definition  of  Constants 


Stage  to  start  computations 

If  0  -  start  computations  at  t0 

If  1  -  start  computations  at  launch 

If  2  -  start  computations  after  motor 
separation 

If  3  -  start  computations  after  airframe 
separation  time 

Stage  to  stop  computations 

If  0  -  stop  computations  at  the  end  of 
the  launcher 

If  1  -  stop  computations  at  the  time  of 
motor  separation 

If  2  -  stop  computations  at  airframe 
separation  time 

If  3  -  stop  computations  in  the  third 
stage  at  altitude  given  in 
exponential  input 

Item  number  of  first  trajectory  to  be 

computed 

Control  for  wind  data 

If  0  -  winds  provided  in  tables  as  f(y) 
will  be  used 

If  1  -  constant  wind  given  in  exponential 
input  will  be  used 

Control  for  stopping  thrust  phase 

If  0  -  stop  on  value  of  t  given  in 
exponential  input 

If  1  -  stop  on  value  of  Vaccel  given  in 
exponential  input 
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Locations 


16,  17,  18 


19,  20,  21 


22,  23,  24 


25,  26,  27 


28,  29,  30 


Definition  of  Constants 

If  2  -  stop  on  value  of  Va  given  in 
exponential  input 

If  3  -  stop  on  value  of  V  given  in 
exponential  input 

Print  control  (output  from  6  degree  of 
freedom  model  used  to  compute  the  thrust 
phase  trajectory) 

If  n  -  (n  any  number  <  1000 )  every  nth 
line  computed  will  be  printed 

Print  control  (output  from  6  degree  of 
freedom  model  used  to  compute  the  after¬ 
burning  phase  of  the  trajectory) 

If  n  -  (n  any  number  <  1000 )  every  nth 
line  computed  will  be  printed 

Print  control  (output  from  particle  model) 
when  used  to  compute  the  afterburning 
phase  of  the  trajectory 

If  n  -  (n  any  number  <  1000 )  every  nth 
line  will  be  printed 

Select  model  for  second  stage  computations 
If  0  -  six  degree  of  freedom  model 
If  1  -  particle  ballistic  model 
Stop  control  (second  stage) 

If  0  -  second  stage  is  terminated  on  input 
value  of  t2E_a,  then  t2E_2  and  enter  a 

determination  of  tg£,  l  ^  ^2E  ^  ^2E— 2 
such  that  third  stage  flight  time  is  minimized 

If  1  -  second  stage  is  terminated  on  input 
value  of  taE 

If  2  -  second  stage  is  terminated  on  input 
value  of  altitude,  ysE 
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Specifications  for  Table  Sizes  and  Formats 

The  coding  permits  flexibility  in  grid  sizes  of  tables  as 
given  below.  Table  sizes  must  be  consistent  with  specifications 
given  on  cards  containing  table  controls.  The  controls  are  as 
follows : 

In  the  format  1414  list 


a.  The  number  of  Mach  numbers  and  angles  of  attack  used  to 
form  the  aerodynamic  tables  Cx(M,a),  Cz0(M,a)>  CZe(M,a), 

c  „„(».<*> 


b.  The  number  of  Mach  numbers  used  in  the  total  given  in 
"a."  above 

c.  The  number  of  angles  of  attack  used  in  the  total  given 
in  "a."  above 


d.  The  number  of  Mach  numbers  and  angles  of  attack  used  to 
form -the  Cmq(M,a)  table 

e.  The  number  of  Mach  numbers  included  in  total  given  in 
"d."  above 


f.  The  number  of  angles  of  attack  included  in  the  total 
given  in  "d."  above 

g.  The  number  of  Mach  numbers  and  angles  of  attack  used  to 
form  the  Cyp(K,a)  table 


h. 


i. 


J- 


k. 


1. 


The  number  of 
"g."  above 

The  number  of 
in  "g."  above 

The  number  of 

The  number  of 
"j."  above 

The  number  of 


Mach  numbers  included  in  the  total  given  in 


angles  of  attack  included  in  the  total  given 

Mach  numbers  and  corresponding  (VH>  values 
Mach  numbers  included  in  the  total  given  in 


Mach  numbers  and  corresponding  C._  (M)  values 

*°A 
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m.  The  number  of  Mach  numbers  included  in  the  total  given  in 
"1."  above 

n.  The  number  of  angles  of  attack  plus  the  total  number  of 
data  points  in  tables  CHOfl,  CNaf>1,  CNOf2,  CN8f2,  Cy7fl , 

cyrf2>  Xyp )  XN1,  and  XN2,  all  assumed  to  be  functions 

of  angle  of  attack  only 

o.  The  number  of  angles  of  attack  included  in  the  total  shown 
in  "n."  above 


p.  The  number  of  Mach  numbers  and  angles  of  attack  used  to 
form  the  aerodynamic  tables  Cx  VZQ)Z>  cz8,2'  Cm,3> 
CmB,2>  ^  cnp>2*  811  assumed  as  functions  of  Mach  number 
and  angle  of  attack 


q.  The  number  of  Mach  numbers  included  in  the  total  given  in 
"p."  above 


r.  The  number  of  angles  of  attack  included  in  the  total  given 
in  "p."  above 

s.  Same  as  "d."  except  that  Cmq(M,a)  applies  to  the  second 
stage  rigid  body  configuration 

t.  Same  as  "e."  except  (^(Mja)  applies  to  the  second  stage 
rigid  body  configuration 


u.  Same  as  "f."  except  Cyp(M,a)  applies  to  the  second  stage 
rigid  body  configuration 

v.  Same  as  "g."  except  Cgp(M)  applies  to  the  second  stage 
rigid  body  configuration 


w.  Same  as  "h."  except  C.a  (M)  applies  to  the  second  stage 

*  8a 

rigid  body  configuration 


x.  The  total  number  of  angles  of  attack  and  corresponding  Cy7 
values 
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y.  The  number  of  Mach  numbers  Included  in  the  total  given 
in  "x."  above 

z.  The  total  number  of  Mach  numbers  and  corresponding 
CX,2(M)  values 

aa.  The  number  of  Mach  numbers  included  in  the  total  shown 
in  "z."  above 

bb.  The  total  number  of  Mach  numbers  and  corresponding 
Cx^3(m)  values 

cc.  The  number  of  Mach  numbers  included  in  the  total  shown  in 
"bb."  above 

dd.  The  number  of  values  of  time  and  corresponding  thrust 

ee.  The  number  of  values  required  to  specify  the  wind  vector 

as  functions  of  altitude  for  the  ascending  and  descending 
branches  of  the  trajectory 

ff.  The  number  of  values  of  time  and  corresponding  nA(t) 

gg.  The  number  of  values  of  time  and  corresponding  p2( t) 

Items  "a."  through  "o."  apply  to  the  first  stage,  "p."  through  "y." 
the  second  stage  rigid  body  trajectory,  "z."  through  "bb."  the 
second  stage  particle  trajectory,  and  "cc."  the  third  stage  particle 
trajectory.  Items  "dd.",  "ff."  and  "gg."  apply  to  the  first  stage 
rigid  body  trajectory  only.  Item  "ee."  applies  to  all  stages. 

The  format  for  all  tables  of  aerodynamic  data  is  8F9.3.  Formats 
for  other  tables  provided  for  are: 

Thrust  8F9 . 3 

Wind  5F14.2 


Air  density  and  2F12.2 

temperature  F12 . 5 

2F12.2 
F12.5 


(time,  thrust,  etc.) 

(altitude,  wind  azimuth,  wind 
magnitude  (ascending)  and  wind 
azimuth,  wind  magnitude 
(descending)) 

(altitude,  air  temperature,  air 
density,  altitude,  air  tempera¬ 
ture,  air  density,  etc.) 
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Input  Generator 

Input  for  as  many  trajectories  as  desired  may  be  generated  from 
a  basic  set  of  input  by  using  the  input  generator.  One  IBM  card  in 
the  format  2413  Indicates  any  one  of  five  types  of  input  generation 
allowed,  the  number  of  changes  to  be  made,  etc.,  as  described  below. 

Let  N1  -  -  -  N8  be  the  eight  words  appearing  on  the  input 
generating  card;  then, 

a.  If  10.  *  1,  ji.e.,  the  first  word  is  001, 

N2  gives  the  number  of  complete  sets  of  input 
to  be  generated 

N3  gives  the  number  of  elements  in  the  exponential 
input  format  to  be  changed 

N4  gives  the  number  of  elements  in  the  integral 
input  format  to  be  changed 

N5  is  000  or  001  depending  on  whether  or  not 
the  input  for  the  basic  trajectory  or  input  from 
the  immediately  preceding  trajectory,  respectively, 
is  to  be  changed 

N6  and  N7  are  zero 

N8  is  000  or  001  depending  on  whether  or  not  compu¬ 
tations  should  cease  or  more  input  follows,  respec¬ 
tively. 

If  N2  is  not  zero  the  card(s)  following  must  con¬ 
tain  words  of  format  4(l3,  E14.7)  which  specify  the  element  number  in 
the  basic  trajectory  to  be  changed  followed  by  the  new  value  for  the 
input  quantity.  If  N3  is  not  zero  the  cards  of  format  4(13,  E14.7) 
are  followed  by  card(s)  of  format  2413  which  contain  the  element 
number  of  the  integral  input  for  the  basic  trajectory  which  is  to  be 
changed  followed  by  the  new  value,  etc. 

b.  If  N1  =  2,  this  signifies  that  a  single  element  is  to  be 

changed  from  the  value  given  in  the  basic  input 

N2  gives  the  element  number  of  the  input  quantity 
in  the  basic  exponential  input  format  to  be 
changed 
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N3  gives  the  number  of  values  the  element  should 
take 

N4  gives  the  element  number  of  the  input  quantity 
in  the  basic  integral  input  format  to  be  changed 

N5  gives  the  number  of  values  the  element  should 
take 

N6  and  N7  are  zero 

N8  is  as  defined  for  the  case  N1  =  1 

If  N2  is  not  zero  the  cards  following  the  input 

generating  card  will  contain  as  many  values  as  are  indicated  by  N3, 

in  the  format  5E14.7.  If  N4  is  not  zero  cards  following  the  input 

generating  card  will  contain  as  many  values  as  are  indicated  by  N5, 

in  the  format  2413. 

c.  If  N1  =  3,  this  signifies  that  a  two  parameter  combination 

of  changes  is  to  be  made  to  the  basic  exponential 
input 

N2  and  N4  give  the  element  numbers  of  the  para¬ 
meters  to  be  changed 

N3  and  N5  give  the  numbers  of  values  assigned  to 
N2  and  N4,  respectively;  values  of  the  parameters 
are  given  in  format  5E14.7 

N5  is  as  defined  for  the  case  N1  *  1 

N6  and  N7  are  zero 

N8  is  as  defined  for  the  case  N1  =  1 

(N2  is  held  constant  while  N3  varies) 

d.  If  N1  =  4,  this  signifies  a  three  parameter  combination  of 

changes  with  N6  and  N7  used  as  described  for  N4 
and  N5  above 

e.  If  N1  ■  5,  this  signifies  that  only  one  trajectory  is  to 

be  computed.  No  input  generation  takes  place 

Item  numbers  automatically  increase  by  one  for  each  succeeding  tra¬ 
jectory  generated. 
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>  If  required 


Arrangement  of  Input  Data  on  Tapes 
Input  Tape  (2) 

Program  and  subroutines 
Exponential  input 
Program  controls 
Wind  table(s) 

Air  density  and  air  temperature  tables 
Special  Tape  (B5)  (First  stage  (six  degree  of  freedom  model)) 
Controls  for  tables 

Mach  number  and  angle  of  attack  grid  for  Cx(M,a),  CZo(M,a), 
CZ8(M,a),  CNb(M,a)  tables 

Cx(M,a)  table 

Czo(M,a)  table 

CZe(M,a)  table 

CNb(M,a)  table 

Mach  number  and  angle  of  attack  grid  for  Cmq(M,a)  table 
Cmq(M>a)  table 

Mach  number  and  angle  of  attack  grid  for  Cyp(M,a)  table 
Cyp(M,a)  table 

Mach  number  grid  followed  by  corresponding  C.  (m)  values 

*p 

Mach  number  grid  followed  by  corresponding  C.R  (M)  values 


a  grid  followed  by  corresponding  values  for  tables  of 

CNOf/a^  CN0f2^a^  CN8f2^a^  Cy7f^a)>  Cy7fs(c*), 

X„b(a),  Xyp(a)’  Xni(a^  311(1  XN2(a) 
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Special  Tape  (B3)  ( Second  stage  ( six  degree  of  freedom  model) ) 

Mach  number  and  angle  of  attack  grid  for  Cx,  Czo,  CZ8,  Cm8, 

C  .  and  C  tables 

np  m 

Cx(M,a)  table 
Cz0(M,a)  table 
CZs(M,a)  table 

Cm8(M>a)  table 
Cnp(M,a)  table 

Cm(M,a)  table 

Mach  number  and  angle  of  attack  grid  for  Cmq(M,Of)  table 
Cmq(M>a)  table 

Mach  number  and  angle  of  attack  grid  for  Cyp  table 
Cyp(M,a)  table 

Mach  number  grid  followed  by  corresponding  C^p(M)  values 
Mach  number  grid  followed  by  corresponding  C^g  (M)  values 

Angle  of  attack  followed  by  corresponding  C  (a)  values 
Special  Tape  (B3)  (Second  stage  (particle  model)) 

Mach  number  grid  followed  by  corresponding  Cx(m)  values 
Special  Tape  (B3)  (Third  stage  (particle  model)) 

Mach  number  grid  followed  by  corresponding  Cx(m)  values 
Thrust  table  (t  followed  by  corresponding  value  of  T,  etc.) 

Wind  table  (y,  Azw^  ascending,  WA  ascending,  Azw^  descending  WA 

descending,  etc.) 


17 


_ NWL  REPCRT  NO.  1812 

Air  temperature  and  air  density  table  (y,  T*  °C,  p,  etc.) 
li  (t)  table  (t  followed  by  corresponding  value  of  ^(t),  etc.) 

li2(t)  table  (t  followed  by  corresponding  value  of  n2(t),  etc.) 
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DISTRIBUTION 

Bureau  of  Naval  Weapons 

DLI-31  4 

RU  1 

RUSD-2  1 

RUSD-22  1 

RUSD-232  1 

Commander 

Armed  Services  Technical  Information  Agency 
Arlington  Hall  Station 
Arlington  12,  Virginia 

Attn :  TIPDR  10 

Commanding  General 

Aberdeen  Proving  Ground 

Aberdeen,  Maryland 

Attn:  Technical  Information  Section 

Development  and  Proof  Services  1 

Attn:  Ballistics  Research  Laboratories  1 

Commander 

Operational  Test  and  Evaluation  Force 

Norfolk  11,  Virginia  2 

Office  of  Technical  Services 
Department  of  Commerce 

Washington  25,  D.  C.  1 

Minneapolis  Honeywell  Regulator  Company 
Duarte,  California 

Via:  INSMAT,  Los  Angeles,  California  1 

Commander 

U.  S.  Naval  Ordnance  Test  Station 
China  Lake,  California 

Attn:  F8053  1 

P8054  1 

4065  1 

50705  1 
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Commanding  Officer 
Picatinny  Arsenal 
Dover,  New  Jersey 
Attn :  ORD-BB 


Librascope  Division,  General  Precision,  Inc. 

Glendale,  California 

Via:  INSMAT,  Los  Angeles,  California 


Commander, 

U.  S.  Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
Attn:  Missiles  Division 


Institute  of  Naval  Studies 
185  Alewife  Brook  Parkway 
Cambridge,  Massachusetts 

Solid  Propellant  Information  Agency 
Applied  Physics  Laboratory 
Johns  Hopkins  University 
Silver  Spring,  Maryland 

Allegany  Ballistics  Laboratory 
P.  0.  Box  210 
Cumberland,  Maryland 

Atlantic  Research  Corporation 
812  North  Fairfax  Street 
Alexandria,  Virginia 

Commander 

U.  S.  Naval  Air  Missile  Test  Center 
Point  Mugu,  California 
Attn:  Technical  Library 

Commanding  General 
Frankford  Arsenal 
Bridge  and  Tacony  Streets 
Philadelphia,  Pennsylvania 

Commanding  General 
Redstone  Arsenal 
Huntsville,  Alabama 
Attn:  Technical  Library 
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DISTRIBUTION  (Continued) 

Commander 

David  Taylor  Model  Basin 

Washington,  D.  C.  1 

Chief  of  Ordnance 
Department  of  the  Army 
Pentagon  Building 
Washington,  D.  C. 

Attn:  ORDTB  -  Ballistics  Section  1 

Director 

National  Bureau  of  Standards 

Washington  25,  D.  C.  1 

General  Electric  Company 

Missile  and  Ordnance  Systems  Department 

Ordnance  Equipment  Operation 

Pittsfield,  Massachusetts  1 

Local: 

K  1 

KB  1 

KBD  5 

KBX  1 

KBX-1  1 

KCPB  1 

ACL  15 

File  1 
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